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We show how clustering as a general hierarchical dynami- 
cal process proceeds via a sequence of inverse cascades to pro- 
duce self-similar scaling, as an intermediate asymptotic, which 
then truncates at the largest spatial scales. We show how this 
model can provide a general explanation for the behavior of 
several models that has been described as "self-organized crit- 
ical," including forest-fire, sandpile, and slider-block models. 
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I. INTRODUCTION 

Clustering and aggregation play an important role in 
many complex systems. In this paper, we present an in- 
verse cascade model for the self-similar growth of clusters. 
Elements are introduced at the smallest scale, which then 
coalesce to form larger and larger clusters. The inverse 
cascade is terminated by the loss of the largest clusters. 
The system is thus in a quasi-steady state with the loss 
of elements in large clusters balanced by the introduc- 
tion of new elements. The clustering process is recog- 
nized to be a branching network similar to a DLA cluster 
or a river network. Individual clusters are analogous to 
branches, and coalescence is equivalent to the joining of 
two branches. 

There is a wide range of applications for this analysis. 
As a specific example, we consider the forest-fire model 
which has been said to exhibit self-organized critical- 
ity M. In one version of the forest-fire model, a square 
grid of sites is considered. At each time step, a model 
tree or a model spark is dropped on a randomly chosen 
site. If the site is unoccupied when a tree is dropped, it is 
"planted." The sparking frequency / is the inverse num- 
ber of attempted tree drops before a spark is dropped. 
If the spark is dropped on an empty site, nothing hap- 
pens; if it is dropped on a tree, it ignites and "burns" 
all adjacent trees in a model forest fire. In this model, 
individual trees are introduced at the smallest scale, clus- 



ters of trees coalesce to form larger and larger clusters. 
Significant numbers of trees are lost only in the largest 
fires that terminate the inverse cascade Q. The noncu- 
mulative frequency-area distribution for the fires is well 
approximated by a power-law relation 



N (X—- 



(1.1) 



with awl. If the sparking frequency / is relatively large, 
the largest fires are relatively small and the self-similar 
inverse cascade is valid only over a relatively small range 
of cluster sizes. If the sparking frequency / is small, 
the fires that terminate the cascade are large and if / 
is sufficiently small the fires will span the entire grid. 
The noncumulative frequency-area distribution of clus- 
ter sizes satisfies equation (1.1) with a « 2 and the cu- 
mulative distribution of clusters with area larger than A 
satisfies equation (1.1) with a « 1. The behavior of the 
one-dimensional forest-fire model has been discussed in 
terms of a cascade by Paczuski and Bak B . The inverse 
cascade analysis is also applicable to the sandpile model 
PI and the slider-block model P . In the sandpile model 
the clusters are the metastable regions that participate 
in avalanches once they are triggered. In the slider-block 
model, the clusters are the metastable regions that par- 
ticipate in slip events once they are initiated. 

One of the most striking patterns in biology is clus- 
ters or aggregations of animals M . Examples range from 
bacteria to whales and include insects, fish, and birds. 
Bonabeau et al. H] showed that the frequency-number 
distribution of whales satisfy equation (1.1) with a k, \. 
The model we present here should also be applicable to 
these biological problems. 



II. HIERARCHICAL CLUSTERING 

We consider a system of stationary entities that we 
shall refer to as elements. In terms of the forest-fire 
model, the elements are the trees that are planted on 



a lattice. The system is growing due to the steady injec- 
tion of new elements that are added to locations that are 
not already occupied by previously injected elements. We 
define connected sets of elements, i.e. groups of elements 
that are in contact, to be clusters. Note, however, that 
our model does not require that elements be confined to 
lattice points. Neighbors can be defined with any metric 
(e.g. distance) condition, or according to a defined graph 
structure (e.g. lattice). In the forest-fire model, clusters 
are the groups of adjacent trees that would burn in a fire 
if a spark dropped on one of the trees in the cluster. We 
construct rules for assigning rank to clusters in such a 
system, based in spirit on the Strahler @] classification 
that was originally developed for branching in river net- 
works. In this classification system, a stream with no up- 
stream tributaries is defined to be of rank one; when two 
rank-one streams combine, they form a stream of rank 
two, and so forth. However, when streams of different 
rank combine, the rank of the dominant stream prevails. 
Our model for the growth of clusters is an extension of 
a scheme developed earlier Q which only allowed for the 
coalescence of clusters of the same rank. The new model 
is much richer in that it accomodates the coalescence of 
clusters of all ranks and can, therefore, describe a much 
wider array of phenomena. 

The rules for our cluster model are: 

1 . We define a single element that is added to a system 
to be a cluster of rank 1. 



5. We terminate the inverse cascade of elements from 
small to large clusters by eliminating clusters of a 
specified high rank. 

In Figure 1, we illustrate how this model works. 

We now wish to establish the dynamical equations gov- 
erning the evolution of this system. Let us define Ni to 
be the number of clusters with rank i, for i > 1. Let rui 
be the average mass — i.e., the number of elements — of a 
cluster of rank i. Then, the total mass Mi of the clusters 
of rank i is given by 



Ali — Nirrij 



(2.1) 



For convenience, we will define the mass of a single el- 
ement to be one, namely mi = 1. For example, in two 
dimensions, we can regard rUi as the mean area Ai of a 
cluster of rank i. This would be the case in the forest-fire 
model. 

We now develop a mean-field approximation describing 
the dynamical evolution prescribed by the mapping rules 
given above. As indicated, we ignore the simultaneous 
coalescence of more than two clusters. We denote the 
instantaneous change in all quantities using the mapping 
symbol i— >. Accordingly, when two clusters of ranks i 
and j coalesce, the values of Ni and Mi are modified as 
follows. For i = j, 



N,, 



iV,+i -1-1, iV, K^ N, 



(2.2) 



2. If a new element is added adjacent to an existing 
cluster, we say that it is added to the cluster with- 
out changing that cluster's rank, unless the cluster 
is a single element. In that special case, we define 
the two elements as forming a cluster of rank 2. 

3. If a new element connects two existing clusters of 
ranks i and j, respectively, then the rank of this 
new cluster is defined as i -I- 1 when i = j and as 
max{i,j} when i ^ j. In words, this is equiva- 
lent to saying that when two clusters of equal rank 
coalesce, then the rank increases by one; however, 
if the two clusters are not of equal rank, then the 
rank of the larger cluster prevails. 

4. If a new element connects three or more clusters, 
then the rank of the new cluster is defined to be 

• the maximal rank of these clusters, when one 
of the clusters has a rank exceeding that of all 
of the others, or 

• the maximal rank of these clusters plus one, 
when there are two or more clusters of the 
same maximal rank. 

[This is a rare event — akin to a four-body 
interaction — and it is ignored in the model equa- 
tions given below.] 



M,+i^Mi+i+2mi, Mi^M,~2m„ (2.3) 

and for i < j, 



N,, 



N, ~ 1, TV, 



TV, 



J' 



Mi 



M, 



Mi i~> Mi — rrii 



(2.4) 



(2.5) 



with equivalent expressions for j < i. In these equations 
for Mj, we have ignored the addition of an element that 
bridges or joins the two clusters. Since mi will be shown 
to increase in an essentially geometric progression with 
respect to the rank i, the omission of that solitary unit 
mass in the calculation does not influence the asymptotic 
properties as i — > oo. 

In our model, coalescence occurs when a new element 
connects two existing clusters. (We have already indi- 
cated that 4-body and higher order effects will be ne- 
glected.) Accordingly, in the mean field approximation, 
we assume that the rate r^ of coalescence between clus- 
ters of ranks i and j is proportional to the product of 
their total numbers, Ni and Nj, and to the product of 
their boundary sizes, £i and ij, and is naturally related to 
the joint probability of the new element connecting two 
pre-existing clusters. For example in two dimensions, 
£i refers to the effective length of the cluster boundary. 
Thus, we assume that 



ry ex NiliNjlj 



(2.6) 



This is an Euclidean approximation, and emerges in the 
spirit of classical kinetic theory, although the mechanics 
of this problem is entirely different. In sections IV and 
VII, this model will be modified to accommodate the 
possible fractal geometry of clusters. 
We now define 



U = Nd, 



(2.7) 



to be the total size of the boundary associated with clus- 
ters of rank i. We select the normalization for our time- 
scale so that rij = LiLj. Accordingly, let C be the in- 
jection rate of single elements, utilizing this time scale. 
The evolution of the system can be determined by appro- 
priately adapting equations (2.2)-(2.5). From equations 
(2.2) and (2.4), we write 



J=2 



(2. 



while the boundary size of individual clusters of rank j 
increase, their absolute numbers will decrease even more 
rapidly so that the total boundary size in clusters of rank 
j will be monotone decreasing. The finite sum, which 
appears in equation (2.11), is somewhat more involved. 
Nevertheless, it is reasonable to expect that the product 
of nik with Lk will steadily diminish as k becomes smaller 
and that negligible contributions emerge from low values 
of k. Finally, it is easy to see that all of the governing rate 
equations will quickly converge, in the sense of an inverse 
cascade from i = 1 to some finite cut-off, as i ^ oo. As 
A^i begins to grow, it provides a stimulus to the growth 
of 7V2j and so on. Similarly, as the masses at each rank 
in the system grow, they will in turn cause the bound- 
ary size Hi of each cluster of rank i to grow, basically in 
proportion to some power in mi. With this intuition in 
hand, we now obtain the steady-state solution for this 
system. 



N,^Li 



2i 



E 



LiLi 



for i > 1. 



(2.9) 



In equation (2.8), we observe that the rate of change in 
the number of clusters of rank 1 is equal to the injection 
rate minus the rate of coalescence of rank 1 clusters to- 
gether with the rate of coalescence of rank 1 clusters with 
clusters of larger rank. The factor of 2 appears because 
two rank 1 clusters were lost in coalescing to form a rank 
2 cluster. Meanwhile, in equation (2.9), we observe that 
the rate of change in the number of clusters of rank i 
is equal to the rate of rank i cluster formation from the 
coalescence of pairs of rank i — 1 clusters, minus the rate 
of coalescence of pairs of rank i clusters, together with 
the rate of coalescence of rank i clusters with clusters of 
larger rank j > i. 

In a similar way, taking into account mi — 1, we can 
express the mass-balance in the system, derived from 
equations (2.3) and (2.5), according to 



Ml = C - 2Ll 



LiLo 



(2.10) 



III. STEADY STATE SOLUTION: CLUSTER AND 
MASS SCALING 

We derive a steady state solution for an inverse cascade 
from equations (2.8) through (2.11). In our inverse cas- 
cade, single elements are introduced at the lowest level, 
and they coalesce to form larger and larger clusters. The 
inverse cascade is terminated by assuming that very large 
clusters are removed from the system. We assume that 
our system develops in a sufficiently large region, so that 
edge effects can be ignored over a long time. Otherwise, 
we will have a completely space-filling solution and per- 
colation effects will govern. We can regard this (limited) 
steady-state solution to be an intermediate asymptotics 
|10{ for our system — our solution will describe the simil- 
itude that emerges before percolation and space-filling 
issues become significant. The steady state solution fol- 
lows when the time derivatives in the left hand sides of 
equations (2.8)-(2.11) vanish with the result 



C = 2LI + J2LiLj. 



(3.1) 



J=2 



M, = 2L^_imi_i 



E 
fc=i 



LiLk-nik - 2Llmi 



Lj_^ =2L^+ J2 ^i^3^ fo'' « > 1- 



3=1+1 



(3.2) 



- ^ L^LjUii, for i > 1. (2.11) 
j=i+i 

Note that equations (2.8) and (2.10) are identical since 
Ml = Ni. 

We observe that the equations above have the poten- 
tial for self-similarity, since most of the sums are infinite 
in extent, and might be expected to be convergent. In- 
tuitively, we expect that Lj will diminish as j increases; 



2L^_^mi_i + y^LiL 



kmk 



k=l 



2Lfmi + ^ LiLjirii, for i > 1. (3.3) 

As noted earlier, equations (2.8) and (2.10) are equiva- 
lent. 



Equation (3.2) has a self-similar solution, since that 
equation is invariant under z i— > i + 1, and depends only 
on Ljl Li. Thus, we seek a solution having the form 



Lj = a.T* 



(3.4) 



where < a; < 1. The first of these constraints on x 
corresponds to boundary sizes being positive, while the 
second is necessary for the summation to exist. We find 
that X satisfies 



2x^ 



J2 X' ' - - = x- 



i+j-2 _ 2j-4 



(3.5) 



j=i+i 



Summing the infinite geometric series explicitly and di- 
viding by a:;^*^'*, we obtain 



2x^ 



l-x 



1, or x'^ -2x^ -x+l = 0. (3.6) 



This equation has a single root in the range < a; < 1, 
namely x ~ 0.55495813... . Given equations (3.1) and 
(3.6), we find that 

C = a^[2 + x/{l- x)] =a^ or a = C^/^ . (3.7) 

Substitution of these results into equation (3.4) gives 

Li = (7^/2 (0.55495813)'^^ . (3.8) 

We now turn our attention to equation (3.3). We sub- 
stitute equation (3.4) into equation (3.3), dividing by 
a?x'^~^ and taking into account equation (3.6). We then 
obtain 

i-l 

2a:;*~^mj_i -|- ^ x^'^'^mk = 

k=l 

™J+2 



2x'+^m^ + 



l-x 



-rrii ~ X rrii 



(3.9) 



This equation does not have an exactly self-similar solu- 
tion, since it is not invariant under i i-^ i + 1. Suppose 
that we make the substitution 



x^ ^rrii — y* ^ 



(3.10) 



assuming that y > 1, whereupon we obtain from sum- 
ming the finite series 



2xy'-^ + a;^ ■ ^' ' ^ =y'-\ 

y- 1 



(3.11) 



We observe that the solution for y in this equation de- 
pends upon i. However, for large i, equation (3.11) 
approximately implies, assuming that we can replace 
y*~^ — 1 by j/*^^, that 2x + x^y/{y — I) — y which we 
rewrite as 



y2 - {x + lfy + 2x = 0. 



(3.12) 



This equation has a unique solution for y > 1, namely 
y — \/x — 1.8019377... . Accordingly, for large i, we have 
asymptotic self-similarity with 



1—7 7 — 1 

rrii ~ oix y 



ac 



,i-l 



(3.13) 



where c = l/x^ = 3.24697602... . With mi == 1, we have 



m, « (3.24697602)' 



(3.14) 



Before moving to issues dealing with fractals and 
branching, the solutions we have just obtained for Li and 
for mi can be immediately exploited. Since Li ex x* and, 
approximately, rrii oc a;^^*, we observe that Li^/rrii ~ 
const. For example in two dimensions, recalling that 
Li = Ni£i and introducing the Euclidean relation that 
ii ex y/rrii, it follows that Nivrii « const, or, equivalently, 
we find the number-mass or number-area relationships 



Ni (X l/rrii ex 1/Ai 



(3.15) 



This is equivalent to equation (1.1) with a — 1. The 
branch numbers Ni are loosely equivalent to a loga- 
rithmic binning of cluster sizes. Logarithmic binning is 
equivalent to a cumulative distribution. Thus, the result 
given in equation (3.15) is in agreement with the distri- 
bution of cluster sizes obtained from the forest-fire model 
as discussed above. The concept of clusters can also be 
extended to both sandpile and slider-block models. In 
these cases, the clusters are the metastable regions that 
will avalanche or slip when an event is triggered. In both 
cases, the cumulative distribution of cluster sizes satisfy 
equation (1.1) with a ss 1. These scaling relationships 
are archetypical of self-organized criticality. Remarkably, 
this scaling has been deduced using solely analytic means 
from our inverse-cascade hierarchical cluster model. 



IV. ADAPTATION FOR FRACTAL PERIMETER: 
CLUSTER AND MASS SCALING 

In the analysis given in the previous sections, we as- 
sumed that the rate of cluster coalescence r^ was pro- 
portional to the linear dimensions of the two clusters as 
given in Equation (2.6). We now generalize this depen- 
dence to account for the possibility of fractal clusters by 
introducing an "efficiency" factor e < 1, with an appro- 
priate scaling such that 



e-^'-^'^NdiNjf.j = e-l^-'lL,Lj 



(4.1) 



As before, rij is the rate of coalescence between clus- 
ters of ranks i and j. This modification can, for exam- 
ple, describe the increased efficiency with which a smaller 
cluster can coalesce with a larger one, since the smaller 
cluster can become attached inside one of the nooks and 
crannies that can characterize a fractal perimeter. 

With this modification, we obtain analogs of equations 
(2.8)-(2.11) 



Ni = Ml = C-2Ll-Y^ e^-^LiLj, 



(4.2) 



j=2 



N,, = Lf 



2L 



E 

J =1+1 



e'^^LiLj, for i > 1, (4.3) 



M, = 2L2_i^".-i + E e*''"'^^^ 



fc^^fc 



fc=i 



2Llmi 



y^ e' ^L^Ljm^, for i > 1. (4.4) 
j=i+i 



In the steady state, we obtain analogs of equations (3.1)- 
(3.3) 



i=2 



(4.5) 



L2_i = 2^2+ ^ e*-JL,Lj, fori>l. (4.6) 

j=i+i 



2xf~^ + x^fi-' ^"^^ 



i-l 



ey-1 



= y 



(4.9) 



We assume that ey > 1. Precisely as in equation (3.11), 
we observe that the solution for y in this equation de- 
pends upon i. However, for large j, equation (4.9) ap- 
proximately implies that 2x + x'^y/{ey — 1) = y which we 
rewrite as 



ey 



(x^ -f 2ex + l)y + 2x = 0. 



(4.10) 



Due to equation (4.8), equation (4.10) has a solution y = 
l/x, for any e. Note that condition ey > I is satisfied for 
y — l/x. Accordingly, for large i, we have asymptotic 
self-similarity with rm « ac*~^, where c = l/x^, as in 
equation (3.13). For example, when e = 3/4, we have 
c = 4. 

It is important to remember that e describes the peri- 
metric fractal scaling for the clusters. The relationship 
between perimetric and areal scaling remains a contro- 
versial topic. However, assuming that one can identify 
an appropriate link between the two, for example in the 
context of forest fire or other models, then the preceding 
discussion makes it possible to identify the frequency- 
area relationship for fractal clusters, in analogy to the 
N oc l/A relationship we identified previously for Eu- 
clidian clusters. 



j-i 



2Ltim,^i+J2e''-'L,L 



krrik 



fe=i 



2L>, -I- ^ e^'^L.Ljm,, for z > 1. (4.7) 

j=2 + l 

Substituting equation (3.4) into (4.6) we obtain an analog 
of (3.6) 



2x^ 



= 1, or x"^ - 2ex^ - x -|- e = 0, 



or e 



1-2x2 



(4.8) 



As the Li are positive, x must be positive from its defini- 
tion in (3.4). Suppose that e — x. Then x — 2x^ — x~x^, 
giving X = 0, a contradiction. Accordingly, for posi- 
tive X, the sign of e — x = x^ / (l — 2x'^) changes only 

at a; = ^1/2 where e changes sign as it passes through 
infinity (due to the denominator). It is easy to see that 
the sign of both e — x and e is positive for < a; < a/1/2, 
and negative for x > \/l/2. As a; < e is nesessary for 
the summation of the geometric series to exist, this im- 
plies that X < •\/l/2. In addition, the condition e < 1 
requires that x < 0.55495813 . . .. For example, x = 0.5 
corresponds to e = 3/4. From equations (4.5) and (4.8), 
we obtain that a = x C^''^ , for any e. 

Let us turn now to the mass balance equation (4.7). 
Substituting (3.4) and assuming x^~^rai = y^~^, we ob- 
tain an analog of equation (3.11) 



V. BRANCHING NUMBERS 

In the analogy between clustering and river networks 
that we have discussed above, we can write for our clus- 
ters 



N,, 



+1 



N, 



(5.1) 



which is known as the bifurcation ratio for river networks. 
Also, we have 



A. 



(5.2) 



which is known as the length-order ratio for river net- 
works. For river networks, the fact that these two ratios 
are almost constant is known as Horton's laws |ll| ]. 

A major step forward in classifying river networks was 
made by Tokunaga |1J]. He extended the Strahler or- 
dering system to include side branching. A first-order 
branch joining another first-order branch is denoted by 
the subscript "11" and the number of such branches is 
iVii, a first-order branch joining a second-order branch 
is subscripted "12" and the number of such branches is 
iVi2; a second-order branch joining a second-order branch 
is subscripted "22" and the number of such branches is 

N22. 

In order to apply the concept of side branching to the 
coalescence of clusters, let us suppose that we have a co- 
alescence of two clusters, of ranks i and j. In the case 



i < j, the cluster of rank i becomes a branch of the cluster 
of rank j. Note that, if the smaller cluster has its own 
branches, these branches are not counted as branches 
of the larger cluster. However, these branches, together 
with all of their branches, etc. are counted as subclus- 
ters of the larger cluster. In analogy to river networks, 
branches are to tributaries as clusters are to drainage 
basins. A branch formed by the cluster of rank i is con- 
sidered to be a subcluster too, and is assigned the rank 
i. Any other subcluster is assigned the rank of a cluster 
from which it first formed as a branch. In analogy to 
river networks, subclusters of a cluster correspond with 
the streams in a drainage basin. The case i > j is treated 
similarly. In the case i ~ j, both clusters of rank i be- 
come branches of rank i of the new cluster of rank i + I. 
Subclusters and their ranks are defined the same way as 
above. 

Let tij be the average number of branches of rank i in 
a cluster of rank j, for i < j, and let n^ be the total 
number of sub-clusters of rank i in a cluster of rank j. 
For i = j, we define tu = na = 1. By definition, for i < j 
we have 



i-i 



= E--^^. 



(5.3) 



k^i 



Moreover, let Ni. 



NjUij be the total number of sub- 



clusters of rank i for all clusters of rank j, and let Ty = 
Njtij be the total number of branches. This classification 
scheme is illustrated in Figure 2. In (a), we have a cluster 
of rank "1" which corresponds to a single tree in the 
forest-fire model. In (b), two clusters of rank "1" have 
coalesced to form a cluster of rank "2." This cluster 
has been joined by a cluster of rank "1." In the forest- 
fire model, two trees on adjacent grid points have been 
joined by a third tree. In (c) and (d), clusters of rank 
"3" and "4" are illustrated. For this example, we have 

«12 = "-23 = «34 = 3, ni3 = n24 == H, "14 = 43, ti2 = 

^23 = hi = 3, ti3 = ^24 = 2, and t^ = 4. 

As before, we regard the coalescence of more than two 
clusters as being exceedingly rare and neglect them in our 
treatment. When two clusters, of ranks i and j coalesce, 
we prescribe the mappings for Nki , Nkj , Tki , and T^ as 
described below. When i = j, 

Nk,i+i'-^ Nk^i+i+2nki, Nki^ Nki-2nki, for fc < i, 

(5.4) 

Ti.i+i 1-^ Tj^i+i + 2, Tk^ii-f Tki-2tki, for /c < i; 

(5.5) 

and when i < j, 

Nkj 1-^ Nkj + Uki, Nki 1-^ Nki - nki, for k <i, (5.6) 



Given the rate of coalescence ry = LiLj, we describe 
the time evolution of the branching process by the fol- 
lowing equations 

Nkj — 2Lj_inkj-i + 2_^ LiLjUki — 2Ljnkj 

i—k 

oo 

- ^ LiLjUkj, for fc < j, (5.8) 
i=j+i 

from equations (5.4) and (5.6), and 

^j-ij = 2I/j_i + Lj^iLj - 2Ljtj^i,j 

oo 

- ^ LkLjtj^i^j, forj>l, (5.9) 

k=3 + l 



J-ij — i^ii^j ^ 2Ljtij — y ^ Lk-Ljtij, tor z < j — 1, 

k=3 + l 

(5.10) 

from equations (5.5) and (5.7). As before, we turn our fo- 
cus to the steady state solution of equations (5.8) through 
(5.10). 



VI. STEADY STATE: BRANCHING NUMBERS 

We begin for the steady state case by setting the time 
derivatives in the left hand sides of equations (5.8)-(5.10) 
to zero. We obtain 



2Lj_inkj-i + 2_^ LiLjUki — 



i^k 



2L'^nkj + ^ LiLjUkj, for k < j. (6.1) 

i=j+i 



2L -_i -I- Lj^iL^ 









1, Tki ^ Tki -tki, ioTk<i. (5.7) 



2L^^tj^i,j+ J2 LkLjt.^i,,, ioTJ>l. (6.2) 
k=j+i 



LiLj — 2L'^tij + y. LkLjtij, for 2 < j — 1. (6.3) 

k=j+i 

We observe that, due to the finite summation present in 
equation (6.1), it is not invariant under j — ki-^j — k+1 
and its solution is not exactly self-similar in j — k. How- 
ever, we now employ the same methodology used in §111 
and obtain asymptotically valid approximate solution. In 
particular, we substitute (3.4) into equation (6.1) and di- 
vide by a'^x^'^'^"^, and we obtain 



i-1 



2x^ ''nfej-i +Y^ 



i-fe+2 
X ^ Uki 



2x^'-'=+2 



=k 
j-k+3 



nkj 



''''kj '^ '^kj • 



(6.4) 



l-x 

Based on our result obtained using equation (3.10), we 
introduce 



Ju I ll^<] ' — ~ /C 



(6.5) 

assuming z > 1, and we obtain from summing the finite 
series in equation (6.4) 

2xz^-^-^ + x^ ■ — — = z^-^. (6.6) 

z — 1 

Approximating z^^^ — \ by z^^^ in the asymptotic limit 
J ^ fc, equation (6.6) approximately implies that 2x + 
x^zjiz —\) — z, or 

z^ - (x + 1)^2 + 2a; = 0. (6.7) 

This latter equation is identical to equation (3.12), and 
has a unique solution z > 1, namely z = Xj x ^ 
1.8019377... and, thereby, demonstrates that the branch- 
ing network description preserves the same structural 
character. Accordingly, for j ^ /c, we have 

nfej « /Jx'^-J'z^-'^ = /3c^-^ (6.8) 

where c = Xjx^ = 3.24697602... as before. Thus, we have 
approximately 

Ukj ~ (3.24697602)^^'' (6.9) 

in the limit j ^ k. For the deterministic example given in 
Fig. 1, we have Ukj ~ 4^"^ for j 3> k. Substituting (3.4) 
into equation (6.3) and dividing by a^a;^^~^, we obtain 

.3 



— ZX Lij ~r 



1 -X 



Uj, for i < j- 1 (6.10) 



which establishes that 



, for I < i — 1. 



(6.11) 



[For the special case that i := j — 1, we have from equa- 
tion (6.2) that 2 + X = tj-i,j-] This, now, is function- 
ally equivalent to the similitude relationship assumed by 
Tokunaga, namely 

tjj = tj_j =aa;'"^'. (6.12) 

Importantly, the behavior that Tokunaga assumed to 
be valid emerges in a completely natural way from the 
underlying mathematics of our inverse cascade. Since 
X = 0.55495813, we have for our inverse cascade 



t. 



(0.55495813) 



i-j+i 



(6.13) 



For the deterministic example given in Fig. 1, we have 



i-j+i 



Uj - (1/2) 

Finally, the connection between our treatment of 
branching and our earlier treatment of clustering needs 
to be established. In particular, we observe that nij turns 
out to be equivalent to ny and that both scale as c-'^^ 
where, as we have already seen, c — 1/x^. 



VII. ADAPTATION FOR FRACTAL 
PERIMETER: BRANCHING NUMBERS 

The branching analysis given in the previous section 
is easily modified to include the fractal perimeter depen- 
dence introduced in equation (4.1). Introducing this re- 
lation into equations (5.8)-(5.10), we obtain 



i-i 



Nkj = 2L]_^nk,j-i + XI <^' ^^»^. 



j ^ki 



i—k 



- 2L^nfcj - X e-' ^LiL^Ukj, ior k < j, (7.1) 

oc 

- Y. ^'^''LkLjt.^i,,, forj>l, (7.2) 



J-ij — £ J^iJ^j ~ 2Ljtij 

oo 

- J2 ^''""LkLjUj, fori<j-l. (7.3) 

k=j+i 

In the steady state, we obtain analogs of equations (6.1)- 
(6.3) 



i-i 



2L]_^nk,j-i+Ye' ^L,L. 



nki 



2L]nkj+ X e^ 'ULjUkj, ior k < j. (7.4) 
i=j+i 



2£j_i +e Lj-iLj - 

OO 

2i2t^_i,, + J2 ^■'"''■^fcijij-ij, forj>l. (7.5) 
k=j+i 



e LiLj — 

oo 

2L]Uj+ Y. ^'^''LkLjUj, ioTKj-1. (7.6) 

k=j+i 

Substituting equation (3.4) into (7.4) and assuming that 
x^~^nkj = z^~^, we obtain, due to (4.8), an analog of 
equation (6.6), namely 



2xz^-''-' + x^"-^ ■ ^'^^' ' ^ - z^ 



ez — 1 



(7.7) 



Assuming ez > 1, we approximate {ezy^^ — 1 by {ezy~^ 
in the asymptotic limit j ^ k. In this case, equation 
(7.7) approximately implies that 2x + x^zj(ez — 1) = z, 
or 



ez^ - [x^ + 2ex + l)z + 2a; = 0. 



(7- 



This latter equation is identical to equation (4.10), and 
has a sohition z = 1/x, for any e. Accordingly, for 
j 3> fc, we have asymptotic self-similarity with Ukj ~ 
I3x^~^ z^~^ = (3c' ~^ , where c = l/x"^ as before. 

Substituting equation (3.4) into equation (7.6) and di- 
viding by a^x^-'"^, we obtain 



txj ^ <i^ ^ i 7 I 



which establishes that 



-U 



tij , for i < j — 1 



U 



^ e«-Ja;»-J+2 



for i < i — \. 



(7.9) 



(7.10) 



[For the special case that i = j — 1, we have from equa- 
tion (7.5) that 2 -I- x/e — ij-ij.] Thus, our modification 
of the Euclidean model to accommodate fractal perimet- 
ric behavior is complete, and the self-similar description 
of the branching process has been shown to follow in a 
completely analogous way. 



VIII. CONCLUSIONS AND DISCUSSION 

In this paper, we have presented an inverse cascade 
model for clustering. This model requires: 

1. The addition of single elements at a prescribed 
small scale; 

2. The consideration of the clustering process as a hi- 
erarchical tree with side branching; 

3. The probability that a cluster of one order will co- 
alesce with another cluster of the same or different 
order is proportional to the product of the number 
of trees of the two orders and the square root of 
their masses (or areas); and 

4. Clusters are lost (destroyed) at a prescribed large 
scale. 

Our inverse cascade model provides a general explana- 
tion for the behavior of several models that have been 
considered to exhibit behavior which has often been de- 
scribed as "self-organized criticality" and occurs in var- 
ious settings including the "forest-fire" model. In this 
model, the planting of individual trees is the introduc- 
tion of single elements, and coalescence occurs when a 
planted tree bridges the gap between two existing clus- 
ters. The model "fires" burn significant numbers of trees 
only in the largest clusters and this terminates the inverse 
cascade. Our model gives the number-mass (or area) dis- 
tribution to be A'' ex \/A] this is also found to be the 
case for the forest-fire model. Our model is also appli- 
cable for the sandpile and slider-block models. In the 
sandpile model, the cluster is the region over which an 



avalanche will spread once it is initiated. In the slider- 
block model, the cluster is the region over which a slip 
event will spread once it is initiated. The initiation of 
an avalanche in the sandpile model and the initiation of 
a slip event in the slider-block model are equivalent to 
a spark being dropped on a tree. In both models the 
clusters grow by coalescence. 

We conclude that these models, which are said to ex- 
hibit self-organized criticality, are neither critical nor self- 
organized. Instead, their behavior is associated with an 
inverse cascade which asymptotically approaches (so long 
as the largest scales are not involved) power-law ("frac- 
tal") scaling. This behavior is related to the self-similar 
direct cascade associated with the inertial-range of fully- 
developed isotropic turbulence. This behavior qualifies 
as a form of "intermediate asymptotics" |10|. It is in- 
teresting to note that earthquakes ||l^, landslides [[l4[ , 
and actual forest fires |15| also have analogous power-law 
frequency- area distributions. 

We have quantified our inverse cascade in terms of a 
branching tree hierarchy with side branching. We have 
adapted the taxonomy used for river networks to the 
growth of our clusters. The order of each cluster is spec- 
ified and, in our mean-field approximation, the number 
of clusters of each order is obtained. We find that this 
distribution is identical to the self-similar side branching 
distribution introduced empirically by Tokunaga. This 
distribution has been found to be applicable for river 
networks flq], DLA clusters |1^, and vein structures of 
leaves iial. 
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FIG. 2. Illustration of the concept of branching applied to 
the coalescence of clusters, (a) A single element, (b) Two 
single elements have been linked to form a cluster or rank 
"2," a third element has joined this cluster as a side branch. 
(c) Two clusters of rank "2" have coalesced to form a cluster of 
rank "3." Another cluster of rank "2" and two single elements 
have been added to this cluster, (d) Two clusters of rank "3" 
have coalesced to form a cluster of rank "4." Another cluster 
of rank "3," two clusters of rank "2" and four single elements 
have been added to this cluster. 
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FIG. 1. Illustration of how two clusters of mass rrii and 
nij coalesce to form a single cluster, when an element (solid 
square) bridges the gap between two clusters. The two clus- 
ters have perimeters £i and £j. This example employs a Carte- 
sian lattice for clarity, although our model does not require a 
lattice structure. 
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